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Combination  of  Satellite  Altimetric  Data  in  the 
Short  Arc  Mode  and  Gravity  Anomaly  Data 


i.  i\ mom  chon 

The  short  arc  adjustment  mode  of  satellite  altimetry  makes  a determination  of 
the  geoid  surface  possible  without  the  requirement  of  highly  precise  reference 
orbits.  In  fact,  approximate  values  of  the  state  vector  parameters  with  some 
fairly  large  a-priori  standard  deviations  are  sufficient.  The  observations  gathered 
over  the  oceanic  regions  represent  the  distances  from  a satellite  to  the  geoid.  For 
many  purposes,  the  geoid  surface  is  assumed  to  coincide  with  the  surface  of 
oceans.  For  a good  determination  of  the  geoid  everywhere,  it  is  necessary  to  com- 
bine satellite  altimetric  data  with  some  other  data  gathered  over  the  continental 
regions,  such  as  gravity  observations.  In  this  paper,  satellite  altimetry  data 
whose  nature  will  be  described  later  is  combined  with  mean  gravity  anomalies  rep- 
resenting over  2200  1°  / 1°  geographic  blocks  covering  many  continental  as  well  as 
oceanic  regions. 

The  pertinent  analysis  of  the  (short  arc)  satellite  altimetry  model  and  the 
gravity  anomaly  model  will  be  preceded  by  a few  explanatory  remarks  with  regard 
to  the  former  model,  the  more  complex  of  the  two.  The  parameters  in  this  model 
are  divided  into  two  basic  groups:  (1)  the  earth  potential  coefficients  which  are 
the  conventional  C's  and  S's  of  the  spherical  harmonic  expansion  including  the 
"central  term"  ro,  and  (2)  six  state  vector  components  for  each  short  arc.  The 

(Received  for  pub  -ation  2 June  1977) 
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parameters  in  the  first  group  may  be  considered  completely  free  to  adjust,  except 
for  the  five  " forbidden"  coefficients,  C'1(),  Cjj,  S.j,  C'^j,  S,^,  which  are  usually 
set  to  zero  and  held  at  that  value.  This  implies  that  the  origin  of  the  Cartesian 
coordinate  system  is  identical  with  the  earth's  center  of  mass  and  that  the  z-axis 
coincides  with  the  earth's  axis  of  rotation.  The  parameters  in  the  second  group, 
together  with  their  a-priori  error  characteristics,  are  assumed  to  be  supplied 
from  an  independent  source.  They  may  be  given  in  one  inertial  coordinate  sys- 
tem for  all  arcs,  or  in  a chosen  inertial  coordinate  system  which  may  vary  from 
arc  to  arc.  For  convenience,  tin-  parameters  may  coincide  with  the  earth  fixed 
(K.  F. ) coordinate  system  at  the  epoch  time,  or  in  the  K.  F.  system  itself.  We 
choose  this  last  possibility  which  lends  to  a whole  adjustment  being  performed  in 
the  E.  F.  system. 

At  this  point,  a discussion  of  what  constitutes  a short  arc  may  be  useful.  As 
pointed  out  in  Ilrown,  * the  rationale  for  the  short  arc  approach  is  that  orbital 
errors  resulting  from  the  enforcement  of  a reasonably  accurate  set  of  potential 
coefficients  truncated  at  fairly  low  degree  and  order  can,  under  proper  circum- 
stances, be  accommodated  by  slight  adjustments  of  the  six  state  vector  parameters. 
For  example,  in  the  case  of  a nearly  circular  orbit  at  an  altitude  of  1000  km,  the 
integration  of  an  arc  of  000  sec  using  enforced  potential  coefficients  (constants) 
truncated  at  (n,  m)  (4,  4)  can  result  in  errors  at  the  extremities  of  the  arc 
approaching  1 to  2 m when  the  state  vector  at  midarc  is  held  to  its  correct  value. 
However,  when  the  state  vector  is  allowed  appropriate  freedom  to  adjust,  the  error 
in  the  resulting  arc  may  be  reduced  to  a level  of  a few  centimeters.  In  the  general 
context  of  the  short  arc  method,  then,  a "short  arc"  is  simply  one  that  is  suffi- 
ciently short  so  that  positional  errors  attributable  to  the  errors  in  enforced  poten- 
tial coefficients  (combined  with  errors  of  truncation)  are  "sufficiently  well"  accom- 
modated over  the  arc  by  an  adjustment  of  the  state  vector.  What  constitutes  "suf- 
ficiently well"  depends  on  the  character  of  the  particular  problem.  By  virtue  of 
such  considerations,  it  can  be  said  that  the  short  arc  method,  by  definition,  does 
not  require  that  specific  attention  be  paid  to  errors  in  potential  coefficients  insofar 
as  each  satellite  arc  is  concerned.  This  is  compatible  with  the  functional  relation- 
ship appearing  in  equation  (43)  of  Brown.  In  particular,  the  coordinates  of  a sat- 
ellite point  on  a short  arc  are  expressed  as  independent  of  the  adjustable  potential 
coefficients. 


1.  Brown,  D.  C.  (19t>7)  llcview  of  Current  Geodetic  Satellite  Programs  and  Recom- 

mendations  for  Future  Programs,  Report  for  NASA  Headquarters,  Contract 
No.  NASW-14(i9. 

2.  Brown,  D.  C\  (1973)  Investigation  of  the  Feasibility  of  a Short  Arc  Reduction 

of  Satellite  Altimetry  for  Determination  of  the  Oceanic  fieoid.  Report  No. 

A FC R I . -T R - 7 3 -0 9 20,  Air  Force  Cambridge  Research  Laboratories  (I.W), 
Bedford,  MA. 
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The  state  vector  components  for  a short  arc  adjustment  often  contain 
standard  errors  which  serve  subsequently  for  weighting  purposes.  Weighting 
must  be  independent  for  each  arc,  Blaha  . Weights  for  orbital  parameters, 
although  weak,  are  required  in  a short  arc  satellite  altimetric  adjustment; 
without  weights  the  determination  of  the  six  state  vector  parameters  would  be 
impossible  since  altimetric  measurements  contain  little  or  no  information 
about  the  spatial  orientation  of  the  orbital  plane.  However,  the  short  arc 
approach  to  the  utilization  of  satellite  altimetry  is  not  dependent  on  intensive 
external  tracking  and  can  be  visualized  as  follows: 

(1)  Reference  orbits  approximately  accurate  to  20  m (rms)  are  initially 
established  from  routine  global  tracking  (VI1F  Doppler  tracking  alone  is  sufficient 
to  establish  suitable  reference  orbits). 

(2)  The  reference  orbits  are  divided  into  a large  number  of  subarcs  situated 
over  oceanic  regions;  such  subarcs  are  limited  in  length  to  at  most  1/4  revolution 
and  average  about  1/G  revolution. 

(3)  Each  subarc  is  treated  as  an  independent  orbit  with  the  epoch  at  midarc 
having  a state  vector  subject  to  a-priori  weighting  consistent  with  the  quality  of  the 
reference  orbit. 

(4)  For  each  subarc,  observation  equations  arising  from  satellite  altimetry- 
are  introduced  (such  equations,  of  course,  introduce  a model  representing  the 
oceanic  geoid). 

(5)  The  adjustment  attempts  the  simultaneous  recovery  of  all  geoidal  param- 
eters along  with  revised  estimates  of  the  state  vectors  of  all  orbital  arcs  (which 
may  number  in  the  thousands  for  dense  global  coverage). 

GEOS-3  altimetric  observations  over  a portion  of  the  North  Atlantic  were  pro- 
vided to  AFGI.  by  NASA;  another  set  of  altimetric  observations  over  adjacent  por- 
tions of  the  Indian  and  South  Pacific  Oceans  was  provided  by  the  Naval  Surface 
Weapons  Center  (NSWC).  Figure  1 depicts  one-hundred  and  twelve  passes  of 
GEOS-3  NASA  data  used  in  the  North  Atlantic  regions,  and  one-hundred  and  eleven 
passes  of  GEOS-3  NSWC  data  used  in  the  Indian  and  Pacific  Oceans  area. 

During  the  initial  preprocessing  of  the  GEOS-3  data,  certain  characteristics 
were  identified  for  editing  criteria.  The  first  editing  level  automatically  examined 
altimetic  measurements  for  gross  errors  that  occurred  from  data  handling  proce- 
dures, e.g. , tape  parity  errors  and  measurement  identification  problems.  The 
second  editing  procedure  was  to  examine  the  altimetic  measurements  for  continuity 


3.  Blaha,  G.  (1073)  The  Combination  of  Gravity  and  Satellite  Altimetry  Data  for 
Determining  the  Geoid  Surface.  Report  of  DBA  Systems,  Inc.;  Af'CRL 
Report  No.  75-0347,  Air  Force  Cambridge  Research  Laboratories,  Hanscom 
AFB,  MA. 
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and  to  eliminate  abrupt  point-to-point  changes.  The  final  editing  is  based  on  a 
three  sigma  criteria  when  compared  to  a polynomial  smoothing  function. 

The  positional  a-priori  standard  errors  of  the  state  vectors  were  assumed  to 
be  24  m,  17  m,  and  8 m in  the  in-track,  cross-track,  and  up-track  directions, 
respectively.  The  velocity  standard  errors  were  similarly  assumed  to  be  0.  02 
m sec,  0.02  m/sec,  and  0.01  m/sec.  The  a-priori  standard  error  was  assumed 
to  be  1.0  m for  all  altimetric  measurements.  The  anomaly  values,  positional 
information,  and  the  a-priori  standard  error  of  each  gravity  anomaly  was  read 
from  a tape  supplied  to  AFGL  by  by  Defense  Mapping  Agency  Aeronautical  Center 
(DMAAC). 

2.  \l).|l  STM  INI'  MODELS 
2.1  Satellite  Mtimetr}  Model 

The  geometry  of  this  model  can  be  symbolized  by  the  vector  equation 
II  - R - r 

as  indicated  in  Figure  2,  where  II,  R,  and  r are  the  magnitudes  of  the  above  vec- 
tors. The  letters  R and  r represent  the  adjusted  quantities  that  fulfill  exactly  the 
prescribed  mathematical  model.  The  approximate  values  of  these  quantities  based 
on  some  initial  values  of  parameters  will  be  denoted  as  R°  and  r°.  The  adjusted 
and  observed  satellite  altimetry  values  will  be  denoted  as  II  and  H , respectively. 
These  notations  will  serve  only  when  forming  the  "discrepancy"  terms  (sometimes 
called  constant  terms)  in  observation  equations;  in  other  instances,  especially 
when  forming  the  partial  derivatives,  the  indices  are  dropped.  Thus  the  state 
vector  components  for  an  arc  are  symbolized  by  X,Y,  Z,  X,  Y,  Z and  the  coor- 
dinates of  S a point  on  this  arc,  bv  X , Y , Z (all  in  the  E.  F.  svstem).  The 

5 S S 

epoch  (of  the  state  vector)  is  denoted  by  t and  the  time  of  the  "event  S"  is  t . 

1 J o s 

The  position  of  S for  a given  t - t is  a function  of  the  state  vector,  written  as 

so 

X X (X,  Y,  Z,  X,  Y,  7.)  , 

s s 

Y Y (X,  Y,  Z,  X.  Y,  Z)  , 

s s 

Z /.  (X,  Y,  /.,  X,  Y,  Z)  . 
s s 

In  agreement  with  the  earlier  discussion,  this  position  is  considered  to  be  indepen- 
dent of  the  adjustable  potential  coefficients  r , C's,  and  S's. 
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Figure  2.  Geometry  of  Satellite  Altimetry  in  the  Farth  Fixed  (Geo- 
centric) Coordinate  System 

in  order  to  express  r as  a function  of  potential  coefficients,  the  values  of  h, 
must  he  known;  they  can  be  computed  from  the  position  of  P which,  in  turn,  can  he 
obtained  from  the  position  of  S and  from  the  measured  value  11^.  However,  the 
direction  of  II  is  unknown  and  will  be  approximated  through  the  use  of  an  ellipsoid. 
In  many  respects,  the  best  approximation  to  the  geoid  is  the  geocentric  mean  earth 
ellipsoid;  the  semimajor  axis  (a)  and  the  eccentricity  (e)  of  this  ellipsoid  offer  the 
means  to  approximate  the  direction  of  II  by  the  geodetic  latitude  and  longitude,  o 
and  X (under  the  present  simplification,  OPS  in  Figure  2 would  lie  in  the  meridian 
plane  whose  longitude  is  A).  The  longitude  is  computed  from 

<SX  VXs  • 

while  the  latitude  may  he  obtained  through  an  iterative  process  as 

tg 6 [Zs  (X^  * Y‘V  "|  <N  ■ h)  [NO  - e“)  • h|  , 

2 2 v 1 2 

N a i\  - e sirT  6)  ~ , 

2 2 1 ’’ 

h (X  . + Y ) ~ cos  4>  - N . 


The  position  of  P is  thus  computed  to  a good  approximation  as  follows: 


F 


t 


X X - cos  0 cos  A , 
P s 

A A II1’  cos  0 sin  A , 
P s 

Z Z - II1’  sin  o ; 

P -s 


this  yields  0 front  the  relation 


tfi" 


7.  (X" 

P P 


V1 


Hoth  0 and  A are  seen  to  depend  on  the  state  vector  parameters  through  toe  position 
of  S, 


(i  f)(X,  V,  Z,  X,  V,  '/.)  , <la) 

A A(X,  V,  Z,  X,  Y,  Z)  . (>>’> 


The  values  of  R and  r can  now  be  expressed  as  functions  of  the  adjustable 
parameters,  namely 


H (X2 

rf 

+ a; 

> •>  i •> 

: f z“>  “ 

■»  s 

(2) 

from  Blaha 

, 2 equation 

(4i>a), 

r 

or 

n 

- 

r r 

o 

1 + 

2 <a 

> (('  eos  mA  • 

Z_»  nm 

S .sin  m>)  P , (sin  n) 

nm  nm 

- 

n 2 

m 0 

J_ 

2 

2 

U) 

3 2 

r r cos 
o 

o (kM) 

(3) 

where  u is  the  angular  velocity  of  earth's  rotation  and  kM  is  the  product  of  the 

gravitational  constant  and  the  earth's  mass;  r , C and  S were  already  labeled 

o nm  nm 

as  r C's  and  S's  while  1*  (sin  0)  are  the  conventional  l.egendre  functions.  In 
o’  nm 

practice,  the  series  expansion  in  (d)  is  truncated  at  some  suitable  n X.  The 
parameter  r corresponds  to  the  radius  of  a fictitious  spherical  earth  having  the 
same  potential  as  the  geoid,  namely 

r kM  1 W 
o o 


1 1 


J 


where  \\  is  the  gravity  potential  of  the  geoid.  To  obtain  a starting  value  of  r^  in 
case  it  is  subject  to  adjustment,  one  could  use  1 , the  potential  of  the  mean  earth 
ellipsoid.  From  (2)  and  ( T ) together  with  equations  (1)  it  follows  that 

R R (X,  V,  /.,  X,  Y , 7.)  , (2') 

r r (r  , C's  S's;X,  V,  Z.  X,  V.  Z)  . (3') 

o 

It  may  he  reiterated  that  although  r is  principally  a function  of  the  potential 
coefficients,  it  becomes  dependent  on  the  orbital  elements  through  the  altimeter 
measurement  which  is  considered  perpendicular  to  the  geoid  or,  to  a good  approx- 
imation, perpendicular  to  the  ellipsoid. 

I pon  using  the  approximate  values  of  parameters  and  the  measured  altimetry 
values,  the  discrepancy  terms  (I.)  may  be  computed  in  one  of  the  two  equivalent 
formulas  below  : 


( ) 
r 


(4a) 


or 


or 


d)  - II 


(4  b) 


We  recall  that  II  is  the  distance  from  a satellite  to  the  geoid  (and  not  to  the  ellip- 
soid) and  that  the  ellipsoid  serves  only  to  approximate  the  direction  along  which 
this  distance  is  measured  in  order  to  compute  d.  The  error  in  d caused  by  this 
approximation  would  amount  to  less  than  4 cm  even  if  the  deflection  of  the  vertical 
were  1'  (for  simplicity,  one  can  take  II  ~ 1000  km).  Considering  the  overall  pre- 
cision of  satellite  altimetry,  such  errors  are  negligible.  It  can  be  shown  that 

d I ri."  (r  ■ II)  II  , 


where 


6 - 0 , 

6 being  the  geocentric  latitude  of  S and  0 representing  the  geocentric  latitude  of  P 
as  in  Figure  2.  In  particular,  b is  computed  from 


For  some  purposes,  the  quantity  d could  be  roughly  evaluated  as 
d x 4,8  m sin^  0 . . . per  1000  km  of  H 
For  the  partial  differentation,  the  model 
H - R - r 

may  be  shown  to  be  entirely  sufficient.  In  the  most  general  form,  the  observation 
equations  are  written  as 

V — 1 d (parameters)  + L 

9 (parameters) 

Upon  considering  ( 2 1 ),  (31)  this  is 


V |3R/a(X,  Y,  Z,  X,  Y,  Z)1 


"dX‘ 

dY 

dZ 

dX 

dl' 

_dZ_ 


- [ar/9(r(),  C's,  S's)] 


dr 

o 

dC's 

d.S's 


- (ar/a(X,  Y,  Z,  X,  Y,  Z)1 


'dX' 

dY 

dZ 

dX 

dY 

.dZ. 


(5) 


where  the  differentiation  of  R with  respect  to  the  state  vector  components  is  done 
via  the  chain  rule  applied  to  matrices  as  follows: 


OR 3R  ^-s) 

9(X,  Y,  Z,  X,  Y,  Z)  3(X  , Y , Z ) 0(X  Y Z,  X,  Y,  Z) 
s s s 


With  the  aid  of  (2),  the  first  matrix  on  the  right  hand  side  is  seen  to  be  composed 
of  the  three  directional  cosines.  The  second  (3  y f!)  matrix,  denoted  Mj,  is  given 
as 


M, 


0(X  , Y , Z ) 3(x  , v , 7.  ) 

S S S a ‘ s s 

8(x  y , zg)  a(x,  y,  z,  x,  y,  ?.) 


3(x,  y,  z,,  x,  y,  z.) 
3(X,  Y,  Z,  X,  Y,  Z) 
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Here,  the  lower  ca.se  letters  are  equivalent  to  their  counterparts  (capital  letters), 
except  that  they  refer  to  a chosen  inertial  system.  The  first  (3  ■ 3)  matrix  on  t ho 
right-hand  side  of  the  last  equation  represents  a rotation  and  the  second  (.1  ■ <:) 
matrix  consists  of  the  first  three  rows  of  the  matrizant.  Upon  using  the  advanta- 
geous choice  of  the  inertial  system  for  each  arc  so  that  it  coincides  with  the  U.  P. 
system  at  the  epoch  (t()),  one  can  express 


9(X.,Y  , Z.) 
s s s 

a(x v , z ) 

b b a 


cos  Mt  - t ) I sin  [u(t  - t ) I 0 
so  So 

- sin  [w  (t , - t ) | cos  |u)(t , - t ) | 0 


"i  o' 

' 0 

-uj 

o' 

T I 

■ T 

u 

0 

0 

. 0 

0 

0. 

and 

9(x,  y,  z,  x,  y,  z) 

9(X,  Y,  Z,  X,  Y,  Z) 

the  partial  derivatives  of  R are  thus  written  as 

9R/9(X,  Y,  Z,  X,  Y,  Z)  [X  / R,  Y /R,  Z /RlM,  . (li) 

s'  s'  ’ s 1 

From  (3),  the  partial  derivatives  9r/9(r  ,C's,  S’s)  can  be  expressed  in  the 

o 1 

following  manner: 


d r / B r - 
1 o 

(i  + p.  + 1 
1 2 

D)  /Q  , 

(7a) 

dr/dC 
' no 

r (a/r)n  I 
o ' 

>n(sin  0)/ Q , 

(7b) 

Dr/9C 

r (a/r)n 

P„„,(sin  0)  cos 

m A 

Q , ) 

(7c) 

nm 

nm 

m > 0 

0r/a  S 

' nm 

r (a/r)n 
o 

P (sin  n)  sin 

nm 

mA 

Q . ) 

(7d> 

where 

Q 1 * (r  / r)(P„  - 3D,' 2) 
o 2 ' 


P 


1 


or 

^ (a/r)n  S(ti) 


•) 


or 

^ n<a  r)n  S(n) 


P 


n 


(sin  (l) 


S(n)  > (C  cos  mX  + S sin  mX)  P 
/.  nm  nm  i 

m 0 

D u2  r cos“  0/C  , C kM/r“ 

If  the  quantity  D is  replaced  by  zero  in  all  of  the  above  formulas,  one  recovers 
equations  (49)  of  Rlaha,  2 where  the  effect  of  the  earth  rotation  was  not  included. 
With  the  present  notations,  equation  (3)  reads 

r r (1  + P.  + I D)  ; 

o 12 

the  term  (1  + Pj  + D)  already  appeared  in  (7a). 

The  differentiation  of  r with  respect  to  the  state  vector  parameters  can  be 
accomplished  along  the  following  lines: 

-Or /9(X,  Y,  Z,  X,  Y,  Z)  =*  -<9r/9^)(9^/90)(90/94)  [9*/9(Xg,  Yg,  Zg)) 

X (3(Xg,  Y Zg)/3(X,  Y,  Z,  X,  Y,  Z)]  , 

where  the  last  (3  x 6)  matrix  had  been  denoted  asMj  and  where  the  ellipsoid 

approximation  to  the  geoid  has  been  used  for  this  purpose  (0r/9X  = 0).  According 
4 

to  Blaha,  this  expression  may  be  developed  to  read 

-0r/9(X,  Y,  Z,  X,  Y,  Z)  » a e2  sin  <j>  cos  <t>  (9<£/9(Xg,  Y , Zg)l  , 
and  finally, 

-9r/9(X,  Y,  Z,  X,  Y,  Z)  =*  (a/R)e2  I-(Xg/R)  sin2  J,  -(Yg /R)  sin2  J, 

(Z  /R)  cos2  ilM1  . (8) 


4.  Blaha,  G.  (1977)  Refinement  of  the  short  arc  satellite  altimetry  adjustment 
model,  paper  published  in  Bulletin  Geodes ique. 
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liquations  ((>)  and  (8)  may  now  be  substituted  into  ( ; •)  so  that  the  form  of 
observation  equations  is  the  following: 


9 2 — 


9 2 — 


V {(X  ,/R)(  1 - (a/R)  e sin  o I , (V.  H)[l-(a  H)  e sin  b]  , 


9 9 — 


(2  , K)l  1 * (a/R)  e“  cos  | } / M 


dX 

dY 

d2 

dX 

dt 

Ld/i  J 


- (3r/&(r  , C's,  S’s)| 


dr 

o 

dCs 

dS’s 


(!>) 


where  L was  given  in  (4a)  or  (4b)  and  where  Or  0(r  Cs,  S’s)  appear  explicitly  in 
(7).  For  each  satellite  arc,  a new  set  of  (weighted)  state  vector  parameters  is 
introduced.  The  adjustment  (of  the  satellite  altimetry  alone  or  in  combination 
with  the  gravity-type  data)  is  effectuated  by  the  method  described  in  Chapter  2 of 
Blaha.  * 

The  values  -(a/R)  e sin-  b,  -(a/R)  e sin-  b,  *(a  R)  e cos  o in  equation  (!') 
may  be  called  "refinement  of  the  partial  derivatives”  as  far  as  the  state  vector 
parameters  are  concerned.  These  corrections  can  change  the  "unrefined"  partial 
derivatives  up  to  0.  (i  percent  of  their  values  (the  average  correction  is  roughly 
0.3  percent).  If  the  corrections  were  not  applied,  the  partial  derivatives  with 
respect  to  the  state  vector  would  often  be  accurate  to  only  2 to  1 significant 
digits.  By  including  them,  the  validity  of  these  partial  derivatives  is  extended 
to  another  2 or  more  significant  digits.  The  remaining  eri  s stemming  from 
the  ellipsoidal  approximation  to  the  geoid  are  believed  to  be  smaller  than  the 
above  corrections  by  at  least  one  order  of  magnitude  and  they  diminish  with 
decreasing  deflections  of  the  vertical.  Furthermore,  such  remaining  errors  are 
random  in  nature  unlike  the  above  corrections.  The  practical  usefulness  of  the 
refinement  in  the  partial  derivatives  may  be  illustrated  as  follows.  \ small  set  of 
satellite  altimetry  data  was  adjusted  by  the  AFGl.  program  SAR11A  (Short  Arc 
Reduction  of  Radar  Altimetry)  described  in  Iladgigeorge  and  Trotter.  ’ In  this 
program,  a different  set  of  coefficients  from  t 's  and  S's  is  used,  but  the  state 
vector  parameters  are  treated  exactly  as  described  earlier.  After  two  iterations 
with  the  "unrefined"  partial  derivatives,  the  parameters  received  no  further 


5.  Iladgigeorge,  Cl.,  and  Trotter,  J.E.  (197ti)  Short  arc  reductions  of  GKoS-3 
altimetry  data.  Paper  presented  at  the  AGU  Fall  Meeting. 


corrections.  However,  using  the  refined  partial  derivatives,  already  the  first 
adjustment  yielded  all  the  final  values.  In  this  way  iterations  may  be  eliminated 
altogether,  which  results  in  important  savings  in  terms  of  computer  run-time. 

2.2  bra'iti  \noimih  Mmlcl 

The  gravity  anomaly  model  used  in  the  present  reductions  corresponds  closely 
to  that  given  by  Rapp,  ' equation  (3). 


Ag  ~ (kM / r“)  (n  - l)(a/r)n  ^ (AC  cos  ml  • AS  sin  mA)  P (sin  d) 

nm  nm  nm 


n 2 


m 0 


(10) 


where  d is  the  geocentric  latitude  of  the  point  associated  with  Ag.  Further, 


AC,0 

C 20 

-C20 

AC40 

C40 

❖ 

-C40 

*C50 

C50 

- C50 

the  other  ACs  and  A.S's  being  essentially  the  conventional  spherical  harmonic 

jjs  sjc  >}e 

potential  coefficients  (C's  and  S's)  themselves.  The  values  C^q,  and  C(;q 

refer  to  the  equipotential  ellipsoid  of  revolution  adopted  as  a base  reference  sur- 
face  in  the  normal  gravity  field.  The  impact  of  replacing  the  coefficients  C^q, 
etc.,  by  zeros  is  negligible.  The  maximum  possible  contribution  of  C(.g  — 
-0.000  000  0051  in  equation  (10)  represents  only  0.03  mgal.  The  contribution 
of  C;jQ  would  be  two  orders  of  magnitude  smaller  and  a similar  pattern  would 
hold  for  further  coefficients.  By  comparison,  the  maximum  contribution  of 
C*Q  ~ +0.  000  0024  is  7 mgal. 

In  considering  rQ  to  be  an  adjustable  parameter  as  in  the  preceding  section, 
the  model  (10)  can  be  rewritten  as 


Ag  — C 


-2<r/r*2)  Ar 
o o 


(n  - l)(a/r)n  AS(n) 


n = 2 


(11) 


5.  Rapp,  R.II.  (1972)  Improved  Models  for  Anomaly  Computations  from  Potential 
Coefficients,  Department  of  Geodetic  Science,  Report  \ro.  181,  The  Ohio 
State  University,  Columbus. 
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-T 


where 


Ar 

o 


r 

o 


s)c 

r 

o 


r"  kM/U 
o o 


l’  being  the  potential  of  the  reference  ellipsoid;  in  practice,  the  initial  value  of 
Ar is  zero  and  the  parameter  r^  is  often  weighted.  The  remaining  new  notation 
is 


n 

AS(n)  = > (AC  „ cos  ml  + AS sin  ml)  P (sin  <t>)  . 

/ . nm  nm  nm 

m 0 

The  series  in-.ide  the  brackets  in  (11)  may  be  expressed  as 
00  oc 

tn  - l)(a/r)n  AS(n)  .=  ^ (n  - l)(a/r)n  S(n)  - (a/r)^  C,,0  P 2 (sin  <b) 
n 2 n- 2 

- 3 (a/r)4  C40  P4  (sin  <t>)  - 5 (a/r)r>  Cfi0  P(.  (sin  J)  . 

The  discrepancy  terms  are  obtained  by  the  usual  approach  as 

L = Ag  - Agb  , (12) 

where  Ag  is  a function  of  the  parameters  to  be  adjusted  (see  equation  11)  and  Ag1’ 
is  the  input  - or  "observed"  - value  of  the  gravity  anomaly.  When  the  gravity 
anomaly  model  is  differentiated  with  respect  to  the  parameters  r , C's  and  S's, 
the  value  of  r may  be  considered  constant.  The  observation  equations  have  the 
following  form: 


V [3Ag/a<ro,  C's,  S's)l 


dr 

o 

dC's 

dS's 


+ L 


(13) 
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where,  from  (11), 


)Ag 

&r 

C) 

-2(r 

A') 

r ~)  , 

o 

(14a) 

1 Ag 

ac 

no 

(n 

- 1 )(a  r)n  1 

’ (sin  6)  , 

( 14b) 

>Ag 

i)C 

nm 

(n 

- l)(a  r)n 

P (sin  b)  cos  m A 

nm  ’ 

1 < 14c) 

m > 0 

1 Ag 

as 

nm 

(n 

- 1 )(a  / r)n 

P (sin  J)  sin  mA 

1 (14d) 

2.11  Error  Propagation 

Either  of  the  two  kinds  of  observation  equations  described  in  the  previous  sec- 
tions may  be  adjusted  separately  if  the  data  is  sufficient.  However,  the  most  use- 
ful results  are  obtained  from  their  combination;  the  reader  is  again  referred  to 
Chapter  2 of  Blaha,  ' for  the  algorithms  that  make  the  combination  and  the  solution 
of  all  the  parameters  possible  (including  the  state  vector  parameters  for  each 

satellite  arc).  After  the  adjustment,  the  variance-covariance  matrix  of  r , Os, 

J o 

and  S's  becomes  available  (denoted  £ r , C's  S's).  This  matrix  is  instrumen- 

o 

tal  in  computing  the  a-posteriori  variances  of  r and  Ag  at  chosen  geographic 

locations.  In  the  first  case,  the  values  of  3 r / 3 ( r , Os,  S's)  arc  utilized  and  in  the 

o 

second  case,  the  values  of  3Ag/0(r^,  C's,  S's)  are  needed.  Both  kinds  of  partial 
derivations  have  been  presented  explicitly  in  equations  (7)  and  (14),  respectively. 
The  formulas  giving  either  variance  are  of  course  unchanged  whether  a separate 
adjustment  or  a combined  adjustment  has  been  performed.  They  are  used  to  gen- 
erate a grid  of  a-posteriori  standard  derivations  in  geoid  undulations  and  in  gravity 
anomalies.  Such  grids  may  serve  to  construct  a contour  map  of  standard  deriva- 
tions; similar  grids  lead  to  contour  maps  of  geoid  undulations  and  gravity  anomal- 
ies themselves. 

The  value  of  r is  closely  related  to  the  geoid  undulation  (N)  as  follows: 

N “ N'  - r - r'  (15) 

where  r'  is  the  radial  distance  to  the  reference  ellipsoid  depicted  in  Figure  2. 

It  can  be  shown  that  for  | n|  as  large  as  120  m,  the  error  in  (15)  is  still  negligible; 
in  particular. 


N'  - N|  < 0.  7 mm  . 


1!> 


Figure  Relation  of  the  Value  of  r With  That  of  the  (Jeoid 
Fndulation  (,\) 


Since  r'  is  constant  (independent  of  the  adjustable  parameters),  the  variance  of  N 
is  the  same  as  the  variance  of  r;  it  is  obtained  from  the  law  of  propagation  of  var- 
iances-covariances,  namely 


a [9r/3(r  , C's,  S's)|2  , 

r o r , t 

o 


[9r/9(r  , C's,  S's)l 


T 


(lfi) 


The  gravity  on  the  geoid  is  by  definition 


g Ag  + -> 


(17) 


and  the  pertinent  variances  similarly  are 


°Ag  it,A«  ',<ro> 


C's,  S's)|  l’r  c,s  s,s|9Ag/8(r  C's,  S's)| 


(18) 


■)  being  the  normal  gravity  (independent  of  the  adjustable  parameters). 


I.  M M K KM  \l  HI-: SI  l.ls 

A combined  global  adjustment  of  satellite  altimetry  data  and  gravity  anomaly 
data  has  been  performed  as  described  in  the  Introduction  of  Section  1.  The 
adjustment  model  that  has  served  in  this  task  is  characterized  by  the  presence 
of  spherical  harmonic  potential  coefficients  up  to  and  including  the  degree  and  order 
(14,  14).  The  option  to  include  a-priori  standard  errors  in  the  coefficients  has  been 
exercised  in  the  following  manner: 
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(—  0.  001 1 C 2 0 1 ) , 

<“0-i|c22l)  . 

(-  0.  1 1 S22 1 ) . 

The  five  "forbidden"  coefficients  have  been  held  fixed  at  the  zero  value.  Due  to 
the  scarcity  of  satellite  altimetry  data,  a separate  altimetry  adjustment  in  the 
(14,  14)  geopotential  model  failed  (the  system  of  normal  equations  was  severely 
ill-conditioned).  The  gravity  anomaly  data  has  been  adjusted  alone  and  in  com- 
bination with  the  altimetry  data.  The  present  gravity  anomaly  data  also  exhibits 
gaps;  furthermore,  in  several  areas  the  a-priori  standard  deviation  is  extremely 
large.  This  is  locally  reflected  by  relatively  large  a-posteriori  standard  errors 
in  geoid  undulations  and  gravity  anomalies  in  the  separate  gravity  anomaly  adjust- 
ment as  well  as  in  the  combined  adjustment. 

The  inclusion  of  satellite  altimetric  data  has  significantly  improved  the  external 
reliability  of  the  gravity  anomaly  adjustment.  This  is  due  not  only  to  the  addition 
of  data  in  areas  where  it  has  been  previously  lacking,  but  also  to  the  fact  that  a 
completely  different  and  independent  source  of  data  has  been  allowed  to  play  its 
role  in  the  determination  of  the  earth's  gravity  field.  As  far  as  the  internal  pre- 
cision is  concerned,  significant  improvements  have  been  noted,  especially  in  the 
areas  of  the  thus  added  data.  This  precision  is  represented  by  a-posteriori  stand- 
ard errors  in  geoid  undulations  and  gravity  anomalies.  As  an  example,  the  stand- 
ard error  in  geoid  undulation  in  the  Indian  Ocean  region  (around  = -P>0°  and 
X 120°)  improved  more  than  10  times  (from  8.  3 m to  0.82  m).  The  internal 
precision  is  of  course  invariably  too  optimistic  with  regard  to  the  "real  world" 
since  it  reflects  the  assumption  that  the  geoid  is  exactly  represented  by  the 
spherical  harmonic  expansion  truncated  at  a given  degree  and  order  (any  further 
coefficients  as  well  as  their  standard  errors  are  considered  to  be  exactly  zero). 
However,  the  a-posteriori  standard  errors  can  be  of  great  assistance  in  determin- 
ing the  goodness  of  fit  accomplished  by  a given  reduction,  in  depicting  areas  of 
insufficient  data,  etc. 

The  combined  adjustment  in  the  (14,  14)  geopotential  model  yielded  the  values 
of  geoid  undulations  (N),  gravity  anomalies  (Ag),  a-posteriori  standard  errors  in 
N(<jjj),  and  a-posteriori  standard  errors  in  Ag  ),  computed  all  in  a 10°  x 10° 
geographic  grid.  These  values  have  been  subsequently  used  for  plotting  contour 
maps.  Figures  4 and  5 show  the  contour  maps  for  N (contour  interval  is  10  m) 
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and  Ag  (contour  interval  is  10  mgal),  respectively.  The  areas  of  insufficient  data 
S on  these  maps  have  been  blanked  out.  The  "cut-off  decision  for  this  purpose  has 

been  facilitated  precisely  by  the  contour  maps  of  standard  errors.  A certain  con- 
tour (in  this  case  the  3 m contour)  has  been  chosen  from  the  plot  associated  with 
beyond  which  the  data  has  been  deemed  insufficient.  The  values  of  in  the 
areas  of  week  geoid  determination  rise  sharply  and  are  three  to  ten  fold  larger 
than  the  values  in  all  the  other  parts  of  the  globe. 

The  geoid  contour  map  has  been  compared  with  the  geoidal  maps  appearing  in 
Ifadgigeorge  and  Trotter.  ‘ In  the  areas  of  comparison,  the  shape  of  the  geoid 

[exhibits  remarkable  similarities.  In  addition,  the  numerical  values  of  geoid  undu- 
lations themselves  agree  very  satisfactorily.  It  is  to  be  noted  that  the  geoid  pre- 
sented in  this  reference  had  in  turn  been  compared,  along  four  different  profiles, 
with  a detailed  gravimetric  geoid  (see  Sections  3.2  and  3.3);  the  results  of  these 
comparisons  were  reported  as  excellent.  The  present  reduction  is  thus  verified 
from  independent  sources. 

In  order  to  provide  more  insight  into  the  interdependence  of  the  degree  and 

order  of  truncation  and  the  detail  in  geoidal  features,  a new  adjustment  has  been 

performed,  this  time  in  conjunction  with  a (13,  13)  geopotential  model.  This  has 

yielded  a new  set  of  values  \,  Ag,  o and  a . The  pertinent  contour  maps  for 

•V  ‘Og 

\ and  Ag  appear  respectively  in  Figures  (i  and  7,  in  which  the  same  areas  have 
been  blanked  out  as  previously.  As  expected,  somewhat  lower  standard  errors 
have  been  noticed  compared  to  their  counterparts  in  the  (14,  14)  geopotential  model. 
From  the  last  two  figures,  a certain  amount  of  smoothing  may  be  observed  if  the 
contours  are  compared  to  these  of  the  two  preceding  figures;  this  fact,  too,  was 
to  be  expected. 

t.  COM  l.l  'ION' 

A simultaneous  adjustment  of  real  satellite  altimetric  and  gravity  anomaly  data 
in  conjunction  with  a (14,  14)  geopotential  model  appears  to  have  been  accomplished 
satisfactorily.  The  geoidal  features  extracted  from  this  reduction  reveal  a good 
agreement  with  those  obtained  from  independent  sources.  A comparison  of  in- 
ternal precision  has  demonstrated  the  beneficial  effect  of  adding  altimetric  data  to 
the  existing  body  of  gravity  anomaly  data.  However,  adding  this  relatively  new  and 
completely  independent  source  to  a more  traditional  type  of  data  in  a combined 
adjustment  is  beneficial  also  for  other  reasons,  such  as  increasing  the  external 
reliability  of  the  results.  An  important  fact  that  has  been  confirmed  during  the 
real  data  reductions  — and  not  only  in  some  computer  simulations  — is  that  highly 
accurate  reference  orbits  are  not  a stringent  requirement  in  the  short  arc  mode 
of  satellite  altimetry. 
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The  present  reduction  in  which  only  a part  of  the  existing  satellite  altimetric 


and  gravity  anomaly  data  has  been  utilized  indicate  that  the  global  geoid  recovery 
with  a standard  error  better  than  1 m could  be  achieved  without  undue  difficulty. 
The  outcome  of  this  study  also  points  to  the  desirability  of  further  refining  the 
program  SAGG  so  that  more  detailed  geoidal  features  in  the  regions  with  abundant 
data  may  be  extracted  without  substantially  increasing  the  computer  core-space 
and  run-time  requirements. 
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